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This paper is devoted to the investigation of the radiation energy losses of an ultra- 
relativistic charged particle channeling along a crystal plane which is periodically bent 
by a transverse acoustic wave. In such a system there are two essential mechanisms 
leading to the photon emission. The first one is the ordinary channeling radiation. This 
radiation is generated as a result of the transverse oscillatory motion of the particle in the 
channel. The second one is the acoustically induced radiation. This radiation is emitted 
because of the periodic bending of the particle's trajectory created by the acoustic wave. 
The general formalism described in our work is applicable for the calculation of the 
total radiative losses accounting for the contributions of both radiation mechanisms. We 
analyze the relative importance of the two mechanisms at various amplitudes and lengths 
of the acoustic wave and the energy of the projectile particle. We establish the ranges of 
projectile particle energies, in which total energy loss is small for the LiH, C, Si, Ge, Fe 
and W crystals. This result is important for the determination of the projectile particle 
energy region, in which acoustically induced radiation of the undulator type and also the 
stimulated photon emission can be effectively generated. The latter effects have been 
described in our previous works. 

1. Introduction 

This paper is devoted to the investigation of the radiation energy losses of an 
ultra-relativistic charged particle channeling along a crystal plane which is period- 
ically bent by a transverse acoustic wave. In such a system, there are two essential 
mechanisms leading to the photon emission. The first one is the ordinary channel- 
ing radiation. This radiation is generated as a result of the transverse oscillatory 
motion of the particle in the channel. This radiation mechanism was suggested in 
Refs. 0j§ for linear crystals and later studied in numerous theoretical and exper- 
imental works (see e.g. Refs. i'Bi). The second one is the acoustically induced 
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radiation (AIR)EHZI. The AIR is emitted because of the periodic bending of the 
particle's trajectory created by the acoustic wave. This mechanism is of particular 
interest, because the AIR has all the features of the undulator radiation, including 
the possibility of the stimulated photon emissioro'Ll. 

Since the mechanism of AIR generation was suggested only recently, let us de- 
scribe it here in more details. This mechanism is illustrated in Fig. ^ 
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Fig. 1. Schematic representation of the initially linear planar channel bent by the transverse acous- 
tic wave. The notations are: d is the channel width, a, A are the AW amplitude and wavelength, 
respectively. 

Under the action of a transverse acoustic wave propagating along the z-direction, 
which defines the center line of an initially straight channel (not plotted in the figure) 
the channel becomes periodically bent. Provided certain conditions are fulfilledlj'D, 
the beam of positrons, which enters the crystal at a small incident angle with respect 
to the curved crystallographic plane, will penetrate through the crystal following the 
bendings of its channel. It results in the transverse oscillations of the beam particles 
while travelling along the z axis. These oscillations become an effective source 
of spontaneous radiation of undulator type due to the constructive interference 
of the ohotons emitted from similar parts of the trajectory. As demonstrated in 
Refs. I3'LI, the number of oscillations can vary in a wide range from a few up to a 
few thousands per cm depending on the the beam energy, the AW amplitude and 
wavelength the type of the crystal and the crystallographic plane. In addition to 
the spontaneous photon emission by the undulator, the scheme presented in Fig. 
leads to a possibility to generate stimulated emission. This is due to the fact, 
that photons, emitted at the points of the maximum curvature of the trajectory, 
travel almost parallel to the beam and, thus, stimulate the photon generation in 
the vicinity of all successive maxima and minima of the trajectory. 

As demonstrated in Refs. B0, the AIR can be well separated from the ordinary 
channeling radiation, if certain conditions on the amplitude and the length of the 
acoustic wave are fulfilled. One of the criteria formulated in Refs. I3'LI for the stable 
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work of the AIR type of undulator is the smaU energy loss of the beam of particles 
penetrating through the crystal. 

In this paper we describe general formalism for the calculation of the total 
radiative energy loss accounting for the contributions of both radiation mechanisms. 
We perform such a calculation for the first time. We analyze the relative importance 
of the two mechanisms at various amplitudes and lengths of the acoustic wave and 
the energy of the projectile particle. We establish the ranges of the projectile particle 
energy, in which total radiative energy loss is negligible for the LiH, C, Si, Ge, Fe 
and W crystals. This result is important for the determination of the projectile 
particle energy region, in which acoustically induced radiation is of the undulator 
type and also the stimulated photon emission can be effectively generated. 

An adequate approach to the problem of the radiation emission by an ultra- 
relativistic particle moving in an external field was developed by Baier and Katkov 
and was called by the authors the "operator quasi-classical method" . The details 
of that formalism can be found in Ref. i. We use this formalism to tackle our 
problem. 

For convenience below we enlist the notations used throughout the paper. 

• e, TO, and 7 = s/mc^ are, respectively, the energy, mass, and relativistic factor 
of a projectile, q is its charge measured in units of the elementary charge e, c 
is the velocity of light. 

• Uu, Xu are the AW amplitude and wave length. The undulator parameter 
equals to Pu = 1 in, where = 27ra„/Au < 1. 

• w is the photon frequency, n is the unit vector in the direction of the emission. 

• L is the crystal thickness, T w L/c is the time of flight of the projectile 
through the crystal. 

• c? is the interplanar spacing. It is assumed that d satisfies the condition d <C 

• U ~ U{p) is the interplanar potential, p ~ [—d/2, -\-d/2] is the distance from 
the midplanc. The quantity Uo stands for the maximum value of U{p). 

• C stands for the factor e/{Rmin<}eU^^^) ~ e d/ (Rmin^qeUo), where Rmin = 
(fc^a,j)~^ is the minimum curvature radius of an acoustically bent channel, 
and fc,i = 2n/Xu, U^^^ « 2Uo/d is the maximum gradient of the interplanar 
field. 

• Oc, Ac are, respectively, the amplitude and wave length characterizing the 
channeling motion. The corresponding undulator parameter reads as pc — 
7^c- The ^-parameter is defined as = SttOc/Ac. More details on these 
parameters are given in the text. 
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• a = e^/hc « 1/137 is the fine structure constant, — /{ra^c^) = 2.818 x 
10^ ^■^ cm is the electron classical radius. 



2. Quasi-classical formalism for the radiative energy loss 

The energy losses, A i?, due to the emission of photons by a charged projectile 
moving in an external field are defined as 

Within the framework of the quasi-classical approach! the distribution of the 
energy radiated by an ultra- relativistic particle (of a spin s = 1/2) in given direction 
n and summed over the polarizations of the photon and the projectile, is given by 
the following expression, which is written up to the terms 7~^: 



The functions Lp{ti,t2) and /(ti,t2) equal to 

ip{ti,t2) = ii - i2 - - n • (ri - r2) , (3) 
c 

fi^^M) = + + (4) 

The notations used are ri_2 = r(ti_2), vi_2 = v(ti^2), with r and v standing for 
projectile's radius vector and velocity, respectively. 

Expression looks almost like the classical formulae, although with quantum 
corrections: 

, e huj 

UJ > LU = — UJ, u~ — , (5) 

£ — rujj e — nui 

which take into account the radiative recoil. 

Expressing the photon frequency and the quantity lu' via the dimensionless 
variable u, fioj = eu/{\ + u), Tioj' = ue, and taking into account the relation 
w^dw = (e/?j) u^/(l + u)^ dw, one obtains the following general expression for the 
relative energy losses: 

AE _ g2 .2 

X Jd^nJ^ ^ dtidt2exp(i^(p{ti,t2)) f{ti,t2). (6) 

Let us transform the functions f{ti, t2) and (p{ti, t2) retaining the terms of orders 
up to and omitting the higher-order terms. 
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To evaluating the factor vi • V2/c^ — 1 from (|^) one makes use of Eqs. (A. 14) 
and ( A.15| ). Then, neglecting the term Ai(ii) Ai(t2), one gets 



vi • V2 ^ _ ]__}_( ^y(*i) "y(^2) 
2 



_1«__L_± ( ilviliZ _ iW£Z ) . (7) 



Substituting this expression into (|j) we obtain 

To transform the phase function 1^9(^1, ^2) from let us first write the trajectory 
Y{t) in the following form: 

Y{t)^e,-{ct + ^z{t))+eyy{t). (9) 

For an arbitrary interplanar potential the function Az(t) can be presented in the 
form (see |l|) 

- I £2 ^2- 

_ V ^ T 

where the parameters and are related to the mean-square values of the trans- 
verse velocities of, respectively, the undulator and the channeling motions 



Az(t) = ~ct 



A.(t) , (10) 
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t U di j ' "2 U dl j ■ ^^^^ 

For the undulator motion (i.e. the motion along the centerline of the acoustically 
bent crystal, yu{t) = o-u sin(27rci/A„)) the first relation from ( |Tl| ) produces the 
result ^„ = 27ra„/A„. Analogously, the parameter may be written as S^c = 
2'naclXci where Oc is the amplitude (mean) of the channeling oscillations < d/2 
and Ac is the period (mean) of the channeling oscillatory motion. To estimate the 
magnitude of one can do the following: ^ adicTc) with Tc ^ \/ m^(F fqeUo 
standing for the period of the channeling oscillations. Hence ^ qeUo/s ^ 1. 

The term Az{t) in the order of magnitude equals to A^it) — 0(^^, £,u^c)- It 
contains only oscillatory terms which satisfy the condition Az (t) = if the averaging 
is carried out over the interval AT > A„/c, r^,. 

Now, to write down the term n- (ri — r2) from (||) let us notice that for an ultra- 
relativistic particle the radiation occurs into a narrow cone with the axis along the 
z-direction. The width of the cone is defined by three parameters, 7~^, £,c, and 
and is equal to 

0„,a. ~ max(7-2, ^ ^ (^2) 

The relations established above allow to write down the following expression for 
(p{ti, 12) which explicitly accounts for all the terms of orders 7"^, and ^„: 

vihM) = ipoihM) + Aip{tiM) , (13) 

VoihM) = KV-i(A,(ti)-A,(t2)) , (14) 
c 

Aip[ti,t2) — —0 —0 C0S(/) , (15) 
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where the following notations are introduced: 

1 



272 



I ^ J. J. 



(16) 



When writing ( p^ we took into account that y{ti) — 2/(^2) ^ 0{^u, £,c)ct which 
results in sin 6* {y{ti) - y{t2)) = 9 {y{ti) - 2/(^2))- 
Using Eqs. (||) and dl^) in one gets 



AS 



8^ \n) 



dti 



{l + uY 



f.eu 

exp \ \— ipo{ti,t2) 



1 + (1+m)2 fVyih) Vy{t2)Y , 2(1 +m) 



(17) 



Due to the relation (|12|) the main contribution to the integral over fin = {&, ' 
comes from the region 9^1. Therefore, one may write 



//•oc />2 
dQn exp (i ^ A^{h,t2)) - J ed9 e'"'^'''/^ y 



-icj' Ay^ cos 0/c 



(18) 



For short the notations u' — ue/fi (see (|5|)) and Ay = y[ti) — y{t2) were used. 
The integrals over <j) and are carried out by using the formulae^u 



/•OO 

i^cos0 ^ 27r Jo(z), / exp(ia2;2) Jo(6x)xdx = ^ 



2a 



exp — 1 



"4a 



where Jo (2:) is the Bessel function of order 0. 
Applying these integrals to (Iq) one gets 



h 27ri 
I = exp 

eu T 



.eu 



2c2t2 



When substituting this result into (17) we first introduce the quantity 



ASi)-^z{t2) 1 (y{ti)-y{t2) 



(19) 



(20) 



(21) 



Hence 



AE 



\h) 



T nT 



"'0 



dti dt2 



{i + uY 



exp (iCu) 



l + (l + ti)2 ^i;,(ti) T;,(t2)V ^ 2{l + u) 



(22) 



Let us demonstrate that the principal contribution to ( |22| ) comes from the region 
I CI < 1. To do this we first evaluate the integral over u and then analyze the result. 
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Making a substitution exp i(u — > i sin (^u in the integrand (which does not affect 
the three-fold integral on the right-hand side) one is left with two basic integrals 



Jl 
h 



V? du 



(1 



u du 



{l + {l + uf) sinC^ 



sin^u , 



[l + uf 

the evaluation of which is elementary but lengthly. The result is 



1 - $ ) /(C) 



Jl ^ - 
J2 = -C + CV(C) + 2C5(C) , 

where the functions /(C) and .g(C) stand for the integrals 



C-$) 3(0 



/(C) 



sin C,u 
l + u 



du. 



3(C) 



cos Qu 
l + u 



du. 



i. 



(23) 
(24) 

(25) 
(26) 

(27) 

Using them 



For ICI > 1 the expansions for /(C) and g{C,) one finds in Ref. 
one gets: 

. _ 6 . _ 12 

Jl ~ J2 ~ 

Large values of C correspond to the t- values |r| :s> h/en"^ (see 
expressions (|8|) taken into account the integrand in ( |22| ) behaves as 
the range of large |C| (or |t|) does not contribute effectively to the integral (|2 

The leading contribution to the energy losses (|2^) comes from the range |C| < 1 
which corresponds to 



(28) 

With the 
Thus, 



ii -U\ = T < 



£ K 



(29) 



3. Energy losses in the case of harmonic interplanar potential 

Instead of evaluating the remaining integrals in (|2^) in the general case we will 
consider the harmonic approximation for the interplanar potential. It allows to 
carry out all the calculations explicitly. The final result can be generalized to the 
case of an arbitrary interplanar potential. 

3.1. Evaluation of the formula for the energy losses. 

The harmonic potential, expressed in terms of the relative distance from the 
midplane, y, can be written in the form 

Uiy) = AUo^. (30) 

The coefficient is chosen to satisfy U{ztd/2) = Uo, with Uo standing for the maxi- 
mum value of the potential. 
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Assuming the following strong inequality (see |^) 



(31) 



which is valid almost in all cases, the dependences Az{t) and y{t) read as (compare 
with (|A.16|) and (|A.19| )) 



y{t) = au sin^J. + Uc sm{flct + c/jq), 



(32) 



A,{t) = --^ <; A„e2 sin2f]„i + Ac^c sin(2f7ci + 20o), (33) 



+4 Ac Cu [cos ((rJc + nu)t + (t)o)+ cos ((r2c - ^u)t + 0o)] 



In these formulas the subscripts "u" and "c" indicate that a quantity is related to 
the undulator motion (the index "u") or to the channeling motion (the index "c"). 

The parameters a„, A^,, ^„ = 27ra„/Ai,, and ri„ — 27rc/Au are explained in . 

The quantities characterizing the channeling motion, which are: the frequency 
of the channeling oscillations flc, the wave length of one oscillation Ac, and the 
parameter ^c = are conveniently expressed through the parameter (see 



Eq. (A.l 



2 lie 



Ac 



ird 



2 /i flc 



(34) 



d /i a 

It is also worth noting that expressions (^) and ( |33| ) embrace both the au ^ d 
and au ^ d regions. In the latter case 5„ = 0, and the formulae produce the result 
for the channeling motion in a linear channel. 

To evaluate the integral ( p^ ) let us first analyze the functions {/S.^{ti) — /S.z{t2)) /ct^ 
{y{ti) — 2/(^2))^ /(ct)^ (see (pi|)) and {vy{ti) — Vy{t2))^ jf? from the integrand. 

It can be easily verified that when substituted into (Az(ii) — ^^{^2)) /ct the 
first term from the right-hand side of ( |3^ ) results in the term of the order the 
second one produces the term ~ and the last one gives the ^ ^c £,u term. To 
illustrate this let us consider the contribution of the first term from (p3f): 



16c TT 



sin 2nut\ — sin 2Vlut2 



ti-to 



£,u ^-u 

Sen 



sin Q„T 



1^. 



COSil„(ti + t2) 



<^. (35) 



Next consideration regarding the difference {Az{ti) — Az{t2)) /ct is that for any 
fixed value of t this function is highly oscillatory because of the factors of the type 
cosr2„(<i + ^2), cosr2c(^i + ^2), and/or sin {{flc ± flu){ti + ^2)/2 + <?!>o). It leads to 
the relation Az(<i) — Az{t2) = for any value of r. The leading value of the integral 
( p2| ) will not be changed if in ( pl| ) one omits these highly oscillatory terms. 

By using similar arguments let substitute the function {yiti) — y{t2))^ / {ctY 
from (|2l|) with its non-oscillatory part. Making use of (|3^ ) one gets 



1 ( y{h)-y{t2) 

2 I CT 



1 (y{h)-y{h) 



4 



sm r]u 



4 



smryc 

Vc 



(36) 
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where r]u,c = ^u,c t/2. 

The main contribution to the integral ( p2[ ) comes from the region of small r (see 
1^), therefore, one may use (sin?7„^c/?7«.c)^ = 1 - ill,c/^- 

Thus, the function C defined in (|2l]) can be substituted with 



C — >xT + aT^, (37) 
where the following short-hand notations are introduced 

.-^ a-±^M±m (38) 

The last function to be transformed is {vy{ti) ~ Vy{t2))^ jc? from the integrand 
in (B3). The non-oscillatory part of it reads as 



^ ^".\ t.- -A- >-'-\ f 

t'. (39) 

Substituting (^7|)-(|39|) into (|22|), and introducing the integration variable r one 
obtains 

Ai; . an^ e . u^du f* dr 



f , u^du r dr r- / 

/ /I , ^4 / — exp lu xr-l-ar 

Jo Jo (1 + u)^ J-t T 
l + (l+u)2 ^2^^2+42^2 2(l + «)- 



.age 
47r ?i 



4 ■ J- 

The limits of the integration over r can be extended to ±cxd. Then one is left 
with two integrals. The first one equals to 

/■°° n ^ 2 d [°° 27ri , 

/ drr exp ViuixT + ar ) — — — / dr cos \u{xt + ar ) = — —r-pr Ai (z). 

J-00 ^ ox Jq [iauy/'^ 

(41) 

Here Ai (z) is the derivative of the Airy's function, and the parameter z = 
xu/{ZauY'^. 

The second integral reads as follows 

/■°° , exp [iw(xr + ar3)l ^. dv . / u^X 
= / dr ^ ^ ^ = 2 i / — sin z w -I- — . (42) 



The derivative of H with respect to z reduces to the Airy's function: dH/dz = 
27riAi(z). Hence 

/•oo 

H = -27ri / dt;Ai(w). (43) 



Here the integration constant is chosen to produce iJ = for z = cxd. 

Substituting (^Tj) and ( ^ ) into (^) and changing the variable of integration 
from M to z and, afterwards, integrating by parts the term containing the integral 
/^°° dw Ai(w) one gets 



Ai; i aq^ £ 2 f"" zdz 



2 c h^^^ Jo PHz) 



1 + P^z) 



Ai'(z) + — Ai(z) 



/32(z) ' ' 3 



(44) 
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where /3(z) = 1 + z"^/^ x- 

The parameter x, the value of which plays the crucial role in defining the mag- 
nitude of the energy losses, is defined as 



?173 

X = — 

e 



1/2 



(45) 



To clarify the meaning of x let us for a moment "switch off" the channeling 
motion by putting = 0. Then (omitting the factor \/2) x ~ ?i7'^C«^«/£ = 
h'j'^Pu^u/s- Here the quantity h'y^pu^u if the frequency of the radiated inten- 



sity maximum, ujmax (in the case pu > 1)0. For lu 3> ujmax the intensity dE/duj 
exponentially decreases. If x < 1 then the intensity reaches its maximum in the 
"physical" domain, i.e. fiojmax < £■ The opposite case x > 1 (and, consequently 
^max > £) corresponds to the situation when a projectile can emit photons of all 
frequencies lying within the range huj = [0, e] so that the spectrum intensity never 
reaches the maximum. In this case the radiative energy losses are dominated by 
the radiation of highly energetic photons, Huj ^ e. 

The analogous arguments can be provided to analyze the case of the channeling 
radiation only, i.e. = 0. 

Hence, the formula (^) is a generalization of the definition of x to the case 
when both motions, the undulator and the channeling ones, exist. This expression 
does not contain the cross-term cx ^„ S.c which could have been originated from 
the interference of the photons emitted due to two different types of motion. The 
absence of the cross-term simply reflects the fact that, generally, the frequencies 
of flu and Qc (and, correspondingly, the frequencies of all possible harmonics) are 
incompatible and, thus, the corresponding electromagnetic waves do not interfere. 

The expression (|4j) can be re-written in the form which is frequently usedJ. 
the theory of energy losses due to the synchrotron and/or undulator radiatiorfl^ 
To do this let us introduce the energy losses calculated in the classical limit. This 
limit corresponds to huj/e <^ 1, which, in turn, means that x ^ 1. Neglecting x in 
the integrand in (H) (it results in putting /3(z) = 1) one gets: 



IJii 



/•oc 




/ zAz 




'o 


12 



Ai'(z). (46) 



Here the second term in the integrand was obtained by integrating by parts the 
term proportional to Ai(z) in (^). 

The integrals are evaluated by using the formula0: 



yielding 



dzz^Ai'(z)--:^^-— r(- + i) r(- + -), (47) 



^ L^X'- (48) 



3 c ?i7^ 
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The functional dependence presented by (^8|) coincides with the well-known 
expression^. 

For a projectile positron (which is, actually, of a prime interest) it is convenient 
to denote the coefficients in (^) in another form by making use of the relations 
e/{hcj) = rrieC^ /{Tic) — ajr,^. Hence 



AS 



- — 

3 re 7 



(49) 



Intserting (|4^) into ( [44[) we obtain the formula for the energy losses within the 
framework of the quasi-classical approach: 



A£; /AS 



where $(x) is given by 



' Ai^(^) + — Ai(z) 



/32(z) 



The functions $(x) and x^^{x) are presented in Fig. 



(50) 



(51) 




1 

Fig. 2. Dependences 'I'(x) and x^^(x) 
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3.2. Comparison of the contributions of undulator and channeling mech- 
anisms to the energy radiative losses. 

From (H) (or (|^)) and (|^-(|l|) it is clear that the magnitude of the energy 
losses for a given projectile energy and a crystal length depends solely on the value 
of the quantity see (p5|). Thus, to compare the relative contributions of the two 
types of radiation it is sufficient to analyze the ratio 

^ _ £,u (52) 
By making use of ( ^ ) and recalling the definition C = e/{R,nin QsU^^^) one gets 

K\ 

In a linear channel (A,j = oo, a,j = 0, Rmin = oo) the range of the Oc values 
corresponding to the channeling motion is ac G [— d/2, d/2]. 

For a channel periodically bent by an acoustic wave the range of the ac values, 
for which the stable channeling motion can occur, is narrower jad < ai™"^'' < d/2. 
Let us first establish the magnitude of ac'""^^ for given projectile energy and AW 
wave length and amplitude. 

Inside the channel the motion of the channeled particle is determined by the 
effective potential 

s 

qeUeffiy) ^ qeU{y) ~ —y, (54) 

where R^^ = -R.~J„ sin(27rz/A„) e [— ^mL'^mLl i^ ^^'^ local curvature radius of 
the channel. Written in terms of the dimensionless variable Y — y/d e [—0.5,0.5] 
and in the case of the harmonic intcrplanar potential the quantity Uef / reads as 

Ueff{Y)^iqeUo (y^ -CY sin^^y (55) 

Schematically, the dependences Ueff{Y) calculated in the points of the AW 
maximum (sin(27rz/An) = +1) and minimum (sin(27r2:/AM) — —1) are presented in 
Fig. ^. The particle will be trapped into the channeling mode only if its energy, 
Ey, associated with the transverse motion is less than AqeUo 

i.e. the minimum of two values AqeUo ± C |y |) |y|_Q 5 which are the heights 

of the asymmetric wells corresponding to the effective potential ( |55| ) in the vicinity 
of the AW minima/maxima. Hence, in an acoustically bent channel the range 
jocl < ai™"^' is determined by following two inequalities (see Fig. ||): 

y2_cy < i(l-C)2, Y^ + CY <^{l^Cf, (56) 



which result in 



|ac|<a(— ' = ^(1-C). (57) 
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Y 



Fig. 3. Effective potential Ui,ff{Y) calculated in the points of the AW maximum, C/(_|_){y) = 
iqeUo (Y^ — CY), and minimum, U(^^){Y) = iqeUo {Y^ + CY), (see (^) versus the dimen- 
sionless distance from the centerline, Y = 2y/d S [—0.5,0.5]. E™"^ = ?7(_){— 0.5) = ?7(^)(0.5) is 
the maximum transverse energy for which channeling in an acoustically bent crystal can occur. 
The quantity Ymax defines the maximum value of the parameter Oc, Ymax = 2a^J"'^^"^ /d. Further 
explanations are given in the text. 



Let us assume that a uniform beam is ideally coUimatcd along the centerline 
when entering the crystal, so that the particles of the beam differ only in the value 
of the initial coordinate yo- Then, the parameter 2ai™°'^''/d = 1 — C defines the 
relative part of the beam particles which are trapped into the channeling mode in 
the acoustically bent crystal. 

Equations (|5^ ) and ( p7|) allow to define the ranges of the parameter Oc for 
which either the undulator or the channeling type of radiation dominate in the 
total energy losses. Namely, for 77 > 1, and, consequently, jad < dC, the undulator 
radiation contributes more to the total energy loss. If 77 > 1 jad > dC then the 
channeling radiation plays the dominant role. If jad = C then both mechanisms 
contribute equally to the radiative losses. The value Oc = C < Cc™"^'' can be reached 
only if C < 1/3. In the opposite case, C > 1/3, the parameter rj is greater than 1 
for all Oc values consistent with jST]), and, hence, the losses due to the undulator 
radiation are higher than those due to the channeling one. 

3.3. Estimation of the magnitude of the parameter x 

With (|5^) taken into account the parameter x from (|4^) becomes: 

X=-^xi™"^^ ^4^2 + 4 |, (58) 
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where (see (M) and (A.18D ) 



xi-^ = ^KA]_,/, = ^^ (59) 



is the maximum value of x due to the channehng radiation in the case of a hnear 
channel, i.e. when Oc™"^'' = d/2. 

The factor 2 ^C^ + al/cP in (||) reaches its maximum value of 2 for C = 1 (to 
get this one substitutes Cc with its maximum value given by ( ^7| ) and, afterwards, 
finds the maximum value with respect to C in the range C = [0,1] within which 
the channeling in an acoustically bent crystal can occur). Hence 

(max) [7^ ^7^ 
Xmax^V^Xc - V 2 — . (60) 

£ d 

Let us estimate the magnitude of the right-hand side of (|60|) for a positron. 



£ — 7 0.511 • 10^ eV. By taking into account the definition (A. 18) one gets after 
some simple algebra 

X^a. =4.3x IQ-S J^- (61) 

Here Vq is the magnitude of qeUo measured in eV and is the intcrplanar spacing 
measured in A. In table 1 the values of Vo and d^ are presented for (110) channels 
in various crystals as indicated. The data for C, Si, Fe, Ge and W were taken from 
Ref. i, the and Vo values for LiH were adopted from Ref. 0. The last column of 
the table corresponds to the values of a positron relativistic factor which produce 
Xmax = 1- Fig. Upresents the dependences Xmaxil) obtained for the (110) channels 
of LiH, C, Si, Fe, Ge and W crystals. 

Table 1. The values of dj^, Vo and 7niin calculated for a positron channelling in (110) planar 
channel in C, Si, Ge, F and W crystals. 



Crystal 


(a1 




7 for 




(eV) 


\max — 1 


C 


1.26 


23 


1.27 X 10" 


Si 


1.92 


23 


1.94 X 10^ 


Ge 


2.00 


40 


1.16 X 10^ 


F 


1.02 


70 


3.39 X 10^ 


W 


1.12 


130 


2.00 X 10"' 



One may conclude, based on the values of 7x„„^=i, and on the curves in Fig. Q 
that in the range 7 < 10^ the parameter x can be chosen to satisfy the condition 

X ^ Xrtiax ^ 1 ■ 

3.4. Realistic calculation of the total energy losses by a positron bunch 
channeling in an acoustically bent crystal. 

Formulae (|49|)-(^) allow to calculate the energy loss for the particular trajectory 
(which is specified by the parameter ac) and for arbitrary crystal length L. A more 
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relativistic factor y 

Fig. 4. Parameter Xmax from ( pl[ ) versus the relativistic factor of a positron channeled in (110) 
channels of various crystals as indicated. 

realistic approach must take into account, firstly, the effect of the decrease in the 
beam volume density with the penetration distance, i.e. the dechanneling effect, 
and, secondly, the distribution of the beam particles in Qc- 

3.4.1. Account J or the dechanneling effect. 

Random scattering of the channeling particle by the electrons and nuclei of 
the crystal leads to a gradual increase of the particle energy associated with the 
transverse oscillations in the channel. As a result, the transverse energy at some 
distance from the entrance point exceeds the depth of the interplanar potential 
well, and the particle leaves the channel. This distance is called the dechanneling 
length _Lrf(7,i?). For a given crystal and channel it depends on a positron energy 
(relativistic factor) and on the curvature radius R 

Therefore, for realistic estimations it is sufficient to assume that the crystal 
length L does not exceed id (7, R) and one may calculate the relative energy loss 
forL-Ld(7,i?). 

It can be demonstrated^ that for a crystal bent with a constant curvature radius 
R the dechanneling length ^^(7, R) satisfies the relation 




(62) 
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where Ld{"f, oo) is the dechannchng length of a positron of the same energy in a 
straight channel {R = oo) and 



is the critical (minimal) radius consistent with the channeling condition in a bent 
crystal, "the centrifugal force < the interplanar force"E3. 

In an acoustically bent crystal the curvature = R^nin sin(27rz/A„) is not 

constant. Therefore, it is natural to consider the mean curvature 1/R which is 
obtained by averaging l/|i?(z)| over the undulator period 

^ = — / fc^a„| sinfc„z|dz = - — . (64) 

Then the dechanneling length in an acoustically bent channel is estimated as follows: 

Ldh,R) = il--Cfjadh), (65) 

TT 

where we introduced the reduced dechanneling length ad (7) = Ld{'j, 00)/^. For a 
given crystal and crystallographic plane this quantity depends weakly on 7. Its ex- 
plicit expression, calculated by using the Lindhard approximation for the potential 
of a planar channel, readal^ 



256 axF 

9^ ~ ln(2e//) - l' 



"(7) = 7T-J — TTTTTTTT^- (66) 



Here qtf — 0.8853Zc ^/^ao and / = 16Z°-^ eV are the Thomas-Fermi atomic radius 
and ionization potential, respectively. Zc is the atomic number of the crystal atoms, 
and flo is the Bohr radius. The dependences ad{"f) for various planar channels can 
be found in Ref. 0. 

Taking into account the quantities introduced above one obtains the following 
expression for the relative energy loss due to the electromagnetic radiation emitted 
in a crystal of the length L = Ld{'j,R): 

^ - 2.3 X 10-^ (^7 j:^ adh) (1 - ^ Cr (c' + I) Hx). (67) 

Here the factor 2.3 x 10~^ (l''^o/d^ originates from 2 /^"i^e (see (^9|)) with 

the parameter x taken in the form given by ( |5^ ) and, in turn, Xc™"^' — Xmax/V^ 
is calculated from (|6l|). 

3.4.2. Energy losses averaged over the parameter a^. 

When a bunch of positrons enters the crystal the particles have various values of 
the initial coordinate yo and of the incidence angle 6q between the momentum of the 
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incident particle and the tangent to the crystal centerline a„ sin(fc„2;). Let us assume 
that the bunch is ideally coUimated so that all the particles have = at the 
entrance. Furthermore, let us assume that the particles are uniformly distributed 
in the yo space. Then, when considering the energy losses by those particles of 
the bunch which are trapped into the channeling mode of motion, it is meaningful 
to carry out the averaging of Ai?/e over the Oc values satisfying (p7[). Thus, the 
average energy losses are defined as follows: 

^ ^ 2.3 X 10-^ ('7 j-^ a,{j) (^1 - I ' G (c, xi""^^^) , (68) 

(max) 



where the function G{C^Xc ) is given by 

G (C, xi""^)) - |- r dx (C^ + x^) <i>(x). (69) 



Here x — Qr /d, $(x) is defined in (|l]) and x = V^Xc"'"''^ VC^ 

+ x^ according to 
(p8[). The upper limit of integration equals X(^ = a^^""^^ jd = (1 — C)I2. 

Although the integration over x can be carried out explicitly, the final result is 
rather cumbersome to be reproduced here. It is simplified in the case x ^ 1 whentj 

*(x) = l-^X + 48x^.. (70) 
16 



Then 
G 

where 



[C,x^r'')) = {c^ + ^) - ^;^(™-)^, +48 (xi""-^)' ^2, (71) 



A. = i|(C^ + .^)^/^[5C^ + 2.x^]+^ln^^i±v|!±M|, (72) 
A, = C' + \G^xl + \xt. (73) 



3.4.3. Results of the numerical calculations of the averaged energy losses. 

As it was mentioned in Sec. there are two types of radiation accompany- 
ing the channeling process of an ultra-relativistic particle in an acoustically bent 
channel: the ordinary channeling radiation and the AIR. Both of them belong 
to an undulator type of radiation. It is known from general theory of a pla- 
nar undulator radiation (see e.g. Ref. □) that its frequency-angular distribution 
dE^ (n) / dwdiln is represented by the sets of characteristic frequencies (harmonics) 
cji^) (if = 1,2,3. . .) each of the width r^-^) = {2/N) {uj'^^^/K), where N is the 
number of the undulator periods. The frequencies of harmonics are defined from 
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the relation hcj'--^^ — 47^riif/(2 + with SI and p standing, respectively, for the 
undulator frequency and parameter. 

In the present paper we analyze mainly the case VLcI^u > 1 (see Eq. (|l])) which 
can be achieved by choosing au ^ d (see Eq. ( A.21D ). Hence, the frequencies of 



harmonics characterizing the AIR and the ch anne ling radiation are well separated 
satisfying the condition wi^-'/wi^^ ^ 1 (Refs. It is the AIR mechanism which 

brings the novelty into the problem. Therefore, let us analyze the stability of the 
undulator AIR radiation towards the decrease in a projectile positron energy due 
to the radiative losses. 

^ The frequencies of harmonics of the AIR radiation one calculates from (see Ref. 

where p„ — 7^u, and 6 is the emission angle with respect to the undulator axis 
(which is the centerline of the initially linear channel, see Fig. |l|). The magnitude 
of 9 satisfies to 9 < 9max = wax{7~^,^„} (Ref. 0). 



From ( |74[ ) one gets the following relation between the total energy losses A i?/e 
and the quantity Awl^'' which is the shift of the Kt\i harmonic frequency from its 
unperturbed value uj\!^^ 

A-i") = . ^ < ^ (75) 

2 + 2t^^7^ + Pu e £ 

The spontaneous AIR radiation formed during the passage of a positron through 
an acoustically bent crystal of total length LdijjR) is stable towards the energy 
loss of the positron provided the shift Acji^-* is smaller than the natural line half- 
width The latter is given by vi^^ /2 = (l/Nu) (ujI'^^ / K), where = 
Ldij, R)/ Xu. Therefore, from j?^) one deduces 



— < (76) 
£ - 2KNu ^ ' 

It was estimated in Ref. i and analyzed in more detail in Ref. Q that the 
realistic range of Nu in an acoustically based undulator is = 10 ... 25 and the 
corresponding number of the harmonics emitted via the AIR mechanism is /sT ~ 1. 
Thus, the stability of the AIR radiation will occur if A£^/e < 0.01. 

Figures || represent the dependences Ai?/e versus 7 calculated, according to 



38|)-(69), for a positron channeling in (110) channels of LiH, C, Si, Fe, Ge and 
W crystals and for several values of the parameter C < 1 which characterizes the 
bending of the channel. 

The chosen crystals are commonly used in experiments devoted to the inves- 
tigation of the channeling phenomena and, in addition, this set includes crystals 
composed of light (LiH, C, Si), intermediate (Fe, Ge) and heavy (W) atoms. 

Figures ^ allow to estimate the range of validity of the condition ( [76|) for Nu = 
10 . . . 25. It is seen from the figures that the inequality ( |76| ) is well-fulfilled for 
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7 < 10^ in the case of LiH crystal, 7 < 5 x 10^ for C, Si, and Ge, 7 < 2 x 10^ for 

Fe and W. 

For higher values of 7 the increase of Ai?/e will result in the shift Awl^-' 
greater than the half-width ri/^'''/2 thus smearing the line over the wider range of 
frequencies. In this case it is meanigful to consider not the harmonic acoustic wave 
but rather the one with varying amplitude and period analogously to how it was 
proposedwhen considering undulator radiation formed in the tampered magnetic 
wigglersllB. 

4. Conclusions. 

In this work we have described the general formalism for the calculation of 
the total radiative energy loss accounting for the contributions of both radiation 
mechanisms, i.e. the acoustically induced radiationE0 and the ordinary channeling 
radiation. Our formalism is based on the quasi-classical approach (see e.g. Ref. 
H). We have analyzed the relative importance of ordinary channeling radiation and 
the AIR to the total radiation energy loss at various amplitudes and lengths of 
the acoustic wave and as a function of the energy of the projectile particle. We 
established the ranges of the projectile particle energy, in which the total radiative 
energy loss is negligible for the LiH, C, Si, Ge, Fe and W crystals. This result is 
important for the determination of the projectile particle energy region, in which 
acoustically induced radiation of the undulator type and also the stimulated photon 
emission can be effectively generated. 

We consider our present research as a milestone for the advanced theoretical 
description of the AIR phenomenon. The goal of further investigation in this field 
is to achieve an accurate quantitative descriptioiLof the undulator radiation and of 
the corresponding laser effect. In our recent workEl we have outlined the phenomena, 
which must be thoroughly considered. 

Here we mention only the problem, which is closely connected to the present 
research. Using the formalism very similar to the one described here, it is interesting 
to calculate the total frequency and angular distribution of of photons emitted due 
to the mechanisms of the ordinary channeling radiation and the AIR. This work is 
in progress at the moment and will become the subject of another publication in 
the near future. 
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relativistic factor y 



Fig. 5. Averaged relative energy losses ( |68[ ) for the crystal length L = L^{'f,R) versus the rel- 
ativistic factor of a positron channeling in (110) channels of various crystals: (a) LiH, (b) C 
(diamond), (c) Si, (d) Fe, (e) Ge, (f) W. The curves corresponds to different values of the 
parameter C = e/{Rm.in Sfif^max) indicated. 
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Fig. 5 {Continued) 
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Fig. 5 {Continued) 
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Appendix A 

5. Particle's motion in an acoustically bent channel 

5.1. Approximations 

In this section we outline the approximations which have been used when con- 
sidering both the motion and the radiation of an uhra-relativistic charged particle 
undergoing planar channeling in a crystal bent by means of a transverse, harmonic, 
plane acoustic wave transmitted along the z-direction which coincides with a crys- 
tallographic direction in the initially linear crystal. It is assumed that the crystal- 
lographic planes are equally spaced and being parallel to the (a;y)-plane. 

The ultra-relativistic particle enters the crystal at z = having only the y— 
and z— velocity components, Vyo, Vzo- It is assumed that Vyo ^ Vzo and Vzo ~ c. 
If one neglects random scattering of the particle by the electrons and nuclei of the 
crystal, then the particle's trajectory lies in the (xy)-plane and is subject to the 
joint action of the interplanar force U'{p) and of the centrifugal force due to the 
crystal bending. 

The necessary (but not sufficient) condition for a projectile to be trapped into 
the channeling mode of motion in a bent crystal is < 0c, where O is the entrance 
angle between the particle's velocity and the channel centerline, and Qc is some 
critical angle (the estimates of Qc in the case of bent channel £3.n be found in Ref. 
til). In a linear crystal Qc coincides with the Lindhard's anglella. 

The strong inequality ^ 0c allows to introduce the continuum approximation 
for the interaction potential U between the charged projectile and lattice atoms 
arranged in atomic planes. In our paper this approximation is used to describe the 
equations of motion for the particle in a bent channel. 

5.1.1. Approximations related to the crystal bending 

The shape of a channel centerline in an acoustically bent crystal is described by 

y{z) = a„ sinfc„z, (A.l) 

with fc„ ~ 27r/A„. 

It was demonstrated in Ref. i that the realistic ranges for the AW amplitude 
and wavelength are 

a„ = 10"^ . . 10"*^ cm, A„ = 10"^...10"^ cm. (A.2) 

Therefore, we introduce an approximation by assuming that the following strong 
inequality is fulfilled: 

C„ = 27r^«l. (A.3) 

All the final formulae written below take into account the contributions of the terms 
proportional to and while the terms the order and higher are omitted 
(except for the cases when the terms are the leading ones). 
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It is easily verified that the length of a ccnterline and the intcrplanar spacing in 
linear [L, d) and in acoustically bent {L\ d!) channels are related as follows 

L' = L{l + 0{eu)). d' = dil + OiC))- (A.4) 
The analogous relationship one finds for the distance p between some inner point 



{y, z) of the bent channel and its centerline (see Fig. A.l) 



p = (y-yor + iz-zor = [y-yoY i + 



{y - yof 



{y~yoY {l+0{eu))- (A.5) 



Therefore, when neglecting the terms of the order and higher one may put L' = L, 
d' ~ d, and disregard the difference between p and y — Qu sin(fc„z). 



Eqs. (A.4)-(A.5) allaw us to introduce the constant field approximation for 
the interplanar potentiaH which assumes that: (1) within any bent channel the 
potential U depends only on the variable p which is the distance of the (y, z) point 
from the channel centerline, (2) the explicit dependence of J7 on p in the acoustically 
bent channel is identical to the U{p) dependence in the linear channel. 



o (y,z) 



(y^Zo) 



Fig. A.l. The coordinates used to describe a projectile position in an acoustically bent channel: 
{y, z) are the Cartesian coordinates of the particles, p = (y — J/o)^ + (-^ ~ -^o)^ ~ {!/ — yo) is the 
distance to the channel centerline. 

In a linear channel, where the coordinates y and p are basically the same, the 
potential depends only on the transverse coordinate y and is a periodic function 
with the period d: U = U{y) = U{y + N d)^ where the index N = . . . , —1, 0, 1, . . . 
enumerates the channels. 

When the channel is bent by the AW then the y coordinate becomes y ~ 
ttu sin(fc„z). The displacement along the z axis is proportional to and, thus, 
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may be disregarded. Hence, the magnitude of U in the hncar channel in the point 
{y, z) corresponds to the the magnitude of the potential in the bent channel cal- 
culated at the point {y + a„ sin(fc„z),z). The dependence oi U on y and z in an 
acoustically bent channel is given by 

U{y,z) = U{y), y ^ y - au sm{kuz) (A. 6) 

The constant field approximation means that the explicit dependence U(jj) in the 
bent channel is equivalent to the dependence U{y) in a linear one. The coordinate 
p, in the bent channel up to the terms ^ fu, is equivalent to y. 

5.1.2. Approximations related to the energy of a projectile 

We assume that for an ultra-relativistic particle, which enters the crystal having 
the energy eq = mc^ 70 (here the index "0" indicates that the quantity is measured 
at the entrance) the following condition is valid while it channels in the crystal 

^ « 1. (A.7) 
Here the quantity Uo stands for the depth of the interplanar potential well. For a 



positron {q = 1) typical values for qeUo are ~ 10 . . . 100 eV, so that ( |A.7|) is well 
fulfilled in the ultra-relativistic case. 



The strong inequality (A.7) justifies the classical description of the particle's 



motion in both linear and bent crystals. 

5.2. The equations of motion 

The Hamiltonian function of a relativistic particle moving in a scalar potential 
U is given by 

H = y/{cp)^ +m^c* + qeU{y), (A.8) 



where p = rwyv is the momentum, and the coordinate y is defined in (A. 6). 

This Hamiltonian does not depend on time. Therefore, the total energy of 
the particle, mjc^ + qeU is conserved. Hence, the relativistic factor satisfies the 
condition 7 = 70 (1 — qeU /sq). If one neglects the term q e U/eq then the relativistic 
factor 

is the integral of motion 7 = 70 = const. 

By introducing the variable y in the equation of motion p = —dH/dv one derives 

y = -— ^ -I- sin(fc„rt), (A. 10) 

m7 dy 

z = 0. (A.ll) 
This system was obtained by omitting the terms (x qe U /eq and oc f", n>2. The 



first term on the right-hand side of ( A. 10 ) represents the acceleration due to the 
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action of the interplanar force. The second one is due to the channel bending and can 
be written in the form ^v^ j R which exphcitly indicates the centrifugal acceleration 
(here, w « c for an ultra-relativistic projectile, and = ku^u sin(/c„z) with z ~ ct 



is the curvature of the centerline (A.l)). This term vanishes in the case of a linear 
channel when a„ = and/or A„ — > oo. 



The equation (A.ll) is readily integrated yielding z = ct. The correction to this 
dependence one finds by using the relation 7 = 70 = const in (A. 9): 



vl{t) 



(A.12) 



Let us now estimate the ratio Vy (t) jc? . There are two typical scales for the velocity 
in the y direction. The first one, v"y^ is related to the motion of the projectile along 
the centerline of the channel. The period of this motion equals k, X^jc^ hence, 
v'y^ ~ 2au/{Xu/c). The second characteristic velocity, is connected with the 
particle oscillations inside the channel due to the action of the interplanar field U . 
The period of this oscillations is estimated as Tc ~ 2tx ^fmF[Jqe\J^' ^ tt m^^ j qeU 



so that 



„(2) 



2d/Tc- Hence 



^1 



< 1. 



(A.13) 



These estimates, combined with the relation Vy{t) — y = y + ccos(fcct) (see (A.l) 



and (|A.6| )) produce 

z{t) =ct + Az{t), 
with Az{t) satisfying the equation 



dAz 
dt 



1 



{y{t) + ^uCCOs{kct)y 



(A.14) 



(A.15) 



For an arbitrary function Uiy) the system ( A.10| )-( |A.15 ) can be easily integrated 
numerically by setting the initial conditions y(0), y(0), Az(0). 

The function y{t) describes the motion of the particle with respect to the cen- 



terline (A.l) of the acoustically bent channel. According to ( |A.6| ), the total y{t) 
dependence is obtained by combining y(t) and the term a„ sin(A:„ci). 

5.3. As a specific case: harmonic approximation for the interplanar 
potential 

The case of the harmonic interplanar potential is of a particular interest because 
it allows an analytical solution of the equations of motion. Substituting the function 



( pO| ) into the right-hand side of (A. 10) one gets the well-known equation for a driven 
pendulum. Its solution y(t) reads 



m 



sin {Qct + (j)o) + 



- 1 



sin 



(A.16) 
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1 /2 

where the quantities il„ — 2itc/\u and il^ — {qeU" / dP m^^ are the frequen- 
cies of, respectively, the undulator motion, i.e. the motion along the centerline of 
the acoustically bent channel, and the channeling motion due to the action of the 
interplanar potential. The amplitude of the channeling oscillations and the pa- 
rameter 00 are defined by the initial conditions of the particle entering the crystal. 
As short-hand notation, a stands for the ratio of the frequencies 



with 



2qeUo 



TT d 



< 1. 



(A.17) 



(A.18) 



We remind that the dependence y{t) defines the deviation of the trajectory from 
the channel's centerline. The total y{t) dependence, according to ( A. 6 ), is obtained 
by combining (A. 16) with the term sin(fc„ci). 

Substituting (A. 16) into (A. 15) and integrating the resulting equation one ob- 
tains the z{t) dependence 



z{t) = ct 



1 - 



272 



4 (^2-1)2 



9 2 



sin(2rici + 2(/io) 



167r (^2 _ 1)^ 



sin 2^ut 



(A.19) 



1 



cos ((fic + + ?f>o) COS ((fie — ^ul^ + ?!'o) 



1 



1 



The analytic solution ( A. 16 ) allows to establish several quantitative statements 
about the conditions which must be satisfied to consider the stable channeling mo- 
tion, i.e. y(t) e [— 0.5d, 0.5d], as well as to make estimates of the relative magnitudes 
of the frequencies fijj and fJc. 

The following conditions must be fulfilled. 

(1) \ac\ < d/2. This condition means that the amplitude of the channeling os- 
cillations due to the action of the interplanar potential must not exceed the 
half- width of the channel. 

It is equivalent to the Lindhard's condition which established the maximal 
entrance angle of the particle with respect to the midplane in a linear channel, 
& < &L — {2qeUo/ sf'^'^ , where 9^ is the Lindhard's critical ang le0. 

This relation reflects the fact that in a linear channel the energy Ey associated 
with the transverse motion is less than the potential barrier Uo- It is easy to 
verify that in the case of a harmonic potential this condition is equivalent to 
lad < d/2. 

(2) The condition for the channeling in a bent channel. Driven oscillations (the 



second term on the right-hand side of ( A. 16 )) must not result in the particle's 
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leaving the channel. Hence, the relation a„/ (cr^ — l) < d/2 must be fulfilled. 
This inequality can be written in the form which clearly exhibits the physical 
condition for the possibility of the channeling process in an acoustically bent 



crystal. Recalling the definition (A. 18) one gets 



1 e t^max -Rn 



< 



1 + d/2au 



< 1. 



(A.20) 



The left-hand side of this relation is the ratio of the maximum centrifugal force 
e/-Rmin = e fc^ flu and the maximum interplanar force ^eC/^iax = ^qeUo/d. 
Channeling in an acoustically bent crystal can occur only if this ration is 
less than 1. For a channfil bent with a constant radius such a condition was 
formulated by Tsyganovll3. 

(3) The relationship between the undulator and the channeling frequencies. From 
( |A.17|) combined with (|A.20[) one obtains the ratio of the frequencies fJc and 

2a„ 



cj' > 1 



> 1, 



(A.21) 



which demonstrates that the frequency of the channeling motion is always 
larger than that of the undulator motion. 
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